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Abstract . This report is concerned with the optimization of 
trajectories for coplanar, aeroassisted orbital transfer (AOT) from 
a high Earth orbit ( HEO ) to a low Earth orbit (LEO) . In particular, 

HEO can be a geosynchronous Earth orbit (GEO) . It is assumed that 
the initial and final orbits are circular, that the gravitational 
field is central and is governed by the inverse square law, and 
that two impulses are employed, one at HEO exit and one at LEO 
entry. During the atmospheric pass, the trajectory is controlled 
via the lift coefficient in such a way that the total characteristic 
velocity is minimized. 

First, an ideal optimal trajectory is determined analytically for 
lift coefficient unbounded. This trajectory is called grazing 
trajectory, because the atmospheric pass is made by flying at 
constant altitude along the edge of the atmosphere until the excess 
velocity is depleted. For the grazing trajectory, the lift 
coefficient varies in such a way that the lift, the centrifugal 
force due to the Earth's curvature, the weight, and the Coriolis 
force due to the Earth's rotation are in static balance. Also, the 
grazing trajectory minimizes the total characteristic velocity and 
simultaneously nearly minimizes the peak values of the altitude 
drop, the dynamic pressure, and the heating rate. 

Next, starting from the grazing trajectory results, a real 
optimal trajectory is determined numerically for lift coefficient 
bounded from both below and above. This trajectory is characterized 
by atmospheric penetration with the smallest possible entry angle, 
followed by flight at the lift coefficient lower bound. 



Ill 


AAR-249 


Consistently with the grazing trajectory behavior, the real optimal 
trajectory minimizes the total characteristic velocity and 
simultaneously nearly minimizes the peak values of the altitude 
drop, the dynamic pressure, and the heating rate. 

Key Words . Flight mechanics, astrodynamics , aeroassisted 
orbital transfer, optimal trajectories, nonlinear two-point 
boundary— value problems, sequential gradient-restoration 
algorithm. 
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= drag coefficient; 

= lift coefficient; 

= drag , N ; 

= dynamic pressure, N/m z ; 

= lift-to-drag ratio modulus; 

= local acceleration of gravity, m/ sec ; 

= acceleration of gravity at h = h a , m/sec ; 

= altitude, m; 

= thickness of the atmosphere, m; 

= heating rate, W/m 2 ; 

= lift, N; 

= mass, kg; 

= radial distance from the center of the Earth 
= radius of the Earth, m; 

= radius of the outer edge of the atmosphere, 

2 

= reference surface area, m ; 

= t/t = dimensionless time; 

= running time, sec; 

= velocity, m/sec; 

= circular velocity at r = r a , m/sec; 

= reference velocity, m/sec; 

= angle of attack, rad; 

= path inclination, rad; 

= Earth's gravitational constant, m 3 /sec 2 ; 

= air density, kg/m 3 ; 

= air density at h = h a , kg/m 3 ; 

= reference air density, kg/m 3 ; 
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t = final time, sec; 

m = Earth's angular velocity, rad/sec; 

AV = characteristic velocity, m/sec. 

Subscripts 

0 = entry into the atmosphere; 

1 = exit from the atmosphere; 

00 = exit from the initial orbit; 

11 = entry into the final orbit. 

Superscripts 

• = derivative with respect to dimensionless time; 

= variable computed in an inertial system. 
Acronyms 

AOT = aeroassisted orbital transfer; 

GEO = geosynchronous Earth orbit; 

HEO = high Earth orbit; 

LEO = low Earth orbit; 

SGRA = sequential gradient-restoration algorithm; 
TPBVP = two-point boundary-value problem. 
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1 . Introduction 

Saving propellant weight and increasing the payload are 
among the most important problems of space transportation. 

Orbital transfer from a high Earth orbit (HEO) to a low Earth 
orbit (LEO) can be made more economic if the aeroassisted orbital 
transfer (AOT ) mode is employed. In particular, HEO can be a 
geosynchronous Earth orbit (GEO) . 

In the AOT mode, use is made of the aerodynamic forces in 
order to achieve the proper amount of velocity depletion during 
the atmospheric pass. Here, the intent is to achieve a specified 
apogee following the atmospheric exit, while minimizing the 
overall propellant consumption and keeping the peak heating rate 
within reasonable bounds during the atmospheric pass. 

Aeroassisted orbital transfer is not only important for 
HEO-to-LEO transfer maneuvers, but may prove to be indispensable 
to future planetary flights. In particular, this statement 
refers to lunar return vehicles. Mars exploration vehicles, and 
Mars return vehicles (Refs. 1-3) . Indeed it is known that, for a 
round-trip Earth-to-Mars mission, the total characteristic 
velocity of the AOT mode is about half that of the all-propulsive 
mode . 

While the AOT prospects are clearly bright, to take proper 
advantage of them it is necessary that guidance and control 
systems be designed with care. In turn, it is imperative that 
optimal trajectories be studied (Refs. 4—8) , since they supply 
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the ideal benchmark that guidance trajectories should strive to 
approach (Refs. 9-14). 

In this report, GEO-to-LEO optimal trajectories are studied 
with the understanding that the atmospheric pass is made with the 
engine shut off and is controlled via the lift coefficient so as 
to deplete excess velocity. The entry and exit values of the path 
inclination are free, and the criterion of optimization is the 
total characteristic velocity, which is being minimized. 

With the above ideas in mind, this paper is organized as 
follows. Section 2 presents the system description. In Section 3, 
several optimization problems are formulated; in Section 4, the 
data used in the numerical experiments on optimal trajectories are 
presented . 

In Sections 5-6, an ideal optimal trajectory is determined 
analytically for lift coefficient unbounded. This trajectory is 
called grazing trajectory, because the atmospheric pass is made 
by flying at constant altitude along the edge of the atmosphere 
until the excess velocity is depleted. For the grazing 
trajectory, the lift coefficient varies in such a way that the 
lift, the centrifugal force due to the Earth's curvature, the 
weight, and the Coriolis force due to the Earth's rotation are in 
static balance. Also, the grazing trajectory minimizes the total 
characteristic velocity and simultaneously nearly minimizes the 
peak values of the altitude drop, the dynamic pressure, and the 
heating rate. 
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In Sections 7-8, starting from the grazing trajectory 
results, a real optimal trajectory is determined numerically for 
lift coefficient bounded from both below and above. This 
trajectory is characterized by atmospheric penetration with the 
smallest possible entry angle, followed by flight at the lift 
coefficient lower bound. Consistently with the grazing trajectory 
behavior, the real optimal trajectory minimizes the total 
characteristic velocity and simultaneously nearly minimizes the 
peak values of the altitude drop, the dynamic pressure, and the 
heating rate. Finally, the conclusions are given in 
Section 9 . 

These results have important implications for the design of 
guidance and control systems. Indeed, it appears possible to 
develop a nominal trajectory having a nearly-constant lift 
coefficient and performance close to that of the optimal 
trajectory; this concept is of particular interest for low 
lift-to-drag ratio spacecraft with a narrow lift range (Refs . 11-12 ) . 
Also, it appears that the properties of minimum characteristic 
velocity, minimum peak altitude drop, minimum peak dynamic 
pressure, and minimum peak heating rate are achieved via the same 
control distribution. Hence, one surmises that AOT guidance and 
control systems can be designed by concentrating on the peak 
altitude drop, because all the other performance indexes are 
directly connected with this quantity. 
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2 . System Description 

We consider coplanar , aeroassisted orbital transfer from a 
high Earth orbit (HEO) to a low Earth orbit (LEO). In particular, 
HEO can be a geosynchronous Earth orbit (GEO) . We employ the 
following assumptions: (i) the initial and final orbits are 

circular and equatorial; (ii) two impulses are employed, one at 
the exit from the initial orbit and one at the entry into the 
final orbit; (iii) the gravitational field is central and is 
governed by the inverse square law. The four key points of the 
maneuver are these: point 00, exit from the initial orbit; point 
0, entry into the atmosphere; point 1, exit from the atmosphere; 
and point 11, entry into the final orbit. See Fig. 1. 

The maneuver starts at point 00 in high Earth orbit with a 
tangential propulsive burn having characteristic velocity AV qo ; 
here, the spacecraft enters an elliptical transfer orbit, 
connecting the points 00 and 0; this elliptical transfer orbit 
is such that its apogee occurs at r 0Q . At point 0, the spacecraft 
enters the atmosphere; after traversing the upper layers 
of the atmosphere, it exits the atmosphere at point 1; during 
the atmospheric pass, the velocity of the spacecraft is depleted, 
due to the aerodynamic drag. At point 1, the spacecraft enters 
an elliptical transfer orbit connecting the points 1 and 11; this 
elliptical transfer orbit is such that its apogee occurs at r^.. 
The maneuver ends at point 11 with a tangential propulsive burn 
having characteristic velocity AV^; here ' ^he spacecraft 
enters the low Earth orbit, in that the magnitude of AV^ 
is such that the desired circularization into LEO is achieved. 
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For the atmospheric portion (h < h ) of the trajectory of 

a 

the AOT spacecraft, we employ an Earth-fixed system; for the 

space portion of the trajectory (h > h ) , we employ an inertial 

a 

system. For h < h , we compute the air density using the US 
— ^ 

Standard Atmosphere, 1976 (Ref. 15); for h > h , we assume that 

a 

the air density is zero. 

2.1. Atmospheric Pass . Because the initial and final 
orbits are equatorial, the atmospheric portion of the trajectory 
of the AOT vehicle is also equatorial if the bank angle is zero. 

We employ the following hypotheses: (i) the atmospheric pass is 

made by flying eastward; (ii) the engine is shut off; hence, the 
AOT vehicle behaves as a particle of constant mass; (iii) Coriolis 
acceleration terms generated by the Earth's rotation are considered, 
while transport acceleration terms are neglected; (iv) the 
spacecraft is controlled via the lift coefficient; (v) under extreme 
hypersonic conditions, the dependence of the aerodynamic coefficients 
on the Mach number and the Reynolds number is disregarded. 

2.2. Differential System . With the above assumptions, and 
upon normalizing the flight time to unity, the equations of motion 
in an Earth-fixed system are given by 

h = x [Vsiny] , ( la ) 

V = t [-D/ m - gsiny] , (lb) 

y = t [L/mV + (V/r - g/V)cosy + 2oj] , (lc) 


with 0 < t < 1. In the above equations. 
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= r e + h. 


= y/r = y/(r e + 


h) 2 . 


( 2 ) 


where y denotes the Earth's gravitational constant. In addition, 
the aerodynamic forces are given by 


D = (1/2)C d pSV 2 , L = (1/2)C l pSV 2 , (3a) 

with p = p (h) . In particular, if a quadratic polar is postulated, 
the relation between the drag coefficient and the lift coefficient 
is given by 


C D = K 0 + K 1 C L + K 2 C L 2 - 

2.3. Control Constraint . To obtain realistic solutions, 
the presence of upper and lower bounds on the lift coefficient 
is necessary. Therefore, the two-sided inequality constraint 


C LL - C L 


- C LU ' 


0 < t < 1, (4) 


must be satisfied everywhere along the interval of integration. 

2.4. Transformation Relations . The following equations 
allow one to pass from quantities computed in an Earth- 
fixed system to quantities computed in an inertial system, and 
viceversa : 


h 


h. 


(5a) 


Vcosy = Vcosy + wr, 


(5b) 


Vsiny = Vsiny. 


(5c) 
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While the above equations are nonlinear, they can be solved 
explicitly to yield inertial quantities in terms of Earth-fixed 
quantities, and viceversa . One obtains the direct relations 


h = h, 


_ 2 2 2 
V = /(V + 2(orVcosy + w r ) , 


(6a) 

(6b) 


y = arctan [Vsiny/ (Vcosy + tor)], 


(6c) 


and the inverse relations 


h = h , 


V = /(V^ - 2(orVcosy + io 2 r 2 ) , 


(7a) 

(7b) 


y = arctan [ Vsiny/ ( Vcosy - tor)]. 


(7c) 


2.5. Boundary Conditions . At the entry into the atmosphere 
(t = 0) and at the exit from the atmosphere (t = 1) , certain 
static and dynamic boundary conditions must be satisfied. Specifically, 
at atmospheric entry, we have 


h 0 = h a ' 

r 00 2<2V * 2 ' V> ‘ 2r 00 r a V * 2 + r a 2 V cos2 ^0 = °' 


(8a) 

(8b) 


where h^ is the thickness of the atmosphere and V* = V^= /Tt/r^) = 
7.832 km/sec is a reference velocity, i.e., the circular velocity 
at r = r . In addition, at atmospheric exit, we have 


h l = h a ' 

r H 2 ( 2v * 2 - V) - 2r ll r a V * 2 + r a 2 V COs2 ^ = °‘ 


(9a) 

(9b) 
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2.6. Summary . The relations governing the atmospheric pass 

include the differential system (l)-(3), the control constraint 
(4), and the boundary conditions (8)-(9), in which the inertial 
quantities are related to the Earth-fixed quantities via the 
transformation relations (6) . In this formulation, the independent 
variable is the time t, 0 <_ t <_ 1 . The dependent variables include 
the state variables h(t), V(t), y(t), the control variable C L (t) , 
and the parameter x. 
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3 . Optimal Control Problems 

Subject to the previous constraints, different optimization 
problems can be formulated, depending on the performance index 
chosen. The resulting optimal control problems are either of the 
Bolza type or the Chebyshev type. 

Problem (Pi ) . Minimum Energy. It is required to minimize the 
energy necessary for orbital transfer. A measure of this energy 
is the total characteristic velocity AV, the sum of the initial 
characteristic velocity AV Q0 , associated with the propulsive burn 
from HEO, and the final characteristic velocity AV^, associated 
with the propulsive burn into LEO. Clearly, 

11 = AV = AV qo + AV n , dOa) 

with 

^00 - / <V r 00 )V * - <V r 00>V os V (10b) 

A ^ll = /<r a /r ll )V * " <r a /r ll ) 'T COSY l' <10C> 

Here, the inertial quantities are related to the Earth-fixed 
quantities via the transformation relations (6) . 

Problem (P2 ) . Minimum Peak Altitude Drop. It is required to 
minimize the peak value of the altitude drop. The performance index 
is given by 

12 = max (h-h), 0 <_ t _< 1 . (11) 

t a 

Here, h = 120 km denotes the altitude corresponding to the outer 

3 . 

edge of the atmosphere. 
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problem (P3 ) . Minimum Peak Dynamic Pressure. It is required 
to minimize the peak value of the dynamic pressure. The performance 
index is given by 

13 = max (DP) , 0 <_ t <_ 1, (12a) 

t 

DP = (1/2) pV 2 . (12b) 

Problem (P4 ) . Minimum Peak Heating Rate. It is required to 
minimize the peak value of the heating rate. The performance index 
is given by 

14 =max(HR) , 0 <_ t <_ 1, (13a) 

t 

HR = C/lp/p*) (V/VJ 3 * 07 . (13b) 

3 

Here, p* = 0.3097E-03 kg/m is a reference density (density at 

the reference altitude h. = h /2 = 60 km), V. = V = 7.832 km/sec 

x a a 

2 

is a reference velocity, and the constant C = 282.3 W/cm represents 
the stagnation point heating rate at p = p* and V = V*, based on 
a nose radius of one foot. 
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4 . Experimental Data 

The following data are used in the numerical experiments on 
optimal trajectories. 

Physical Constants . The physical constants used in the 

computations are as follows: the radius of the Earth is r g = 6378 km; 

the radius of the outer edge of the atmosphere is r = 6498 km; 

a 

the thickness of the atmosphere is h = 120 km; the Earth's 

a 

3 2 

gravitational constant is y = 0.3986E+06 km /sec ; the Earth's 
angular velocity is w = 0.7292E-04 rad/sec. 

Atmospheric Model . The atmospheric model used is the US 
Standard Atmosphere, 1976 (Ref. 15) . In this model, the values of 
the density are tabulated at discrete altitudes. For intermediate 
altitudes, the density is computed by assuming an exponential 
fit for the function p (h) . 

Transfer Maneuver . A GEO-to-LEO transfer maneuver is considered. 

The GEO radius is r Q0 = 42164 km, corresponding to r^/r^ = 6.4888. 

The LEO radius is r^ = 6708 km, corresponding to r 11 / r a = 1.0323. 

Spacecraft SI . See Refs. 13-14. The mass per unit reference 

surface area is m/S = 300 kg/m 2 . The drag polar C D (C L ) is shown 

in Fig. 2A and is represented by Eq. (3b), with Kq = 0.10, 

= 0.00, K 2 = 1.11. The lift coefficient is subject to Ineq. (4), 

with C _ = -0.90 and C = +0.90. The maximum lift-to-drag ratio 
LL IjU 

modulus is E - 1.50 and occurs for C T = 0.30. 
max Xj 

In the numerical experiments, the possibility of a smaller 
lift range is considered via three alternative values for the 
lift coefficient lower bound: 


\ 
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(i) C TT = -0.90, (ii) C _ _ = -0.70, (iii) C Tr = -0.50. (14) 

1j.Li LL LL 

For these cases, the maximum lift-to-drag ratio modulus is E max = 1*50 

and occurs for C = 0.30. 

Li 

Spacecraft S2 . See Ref. 11. The mass per unit reference 

2 

surface area is m/S = 117.3 kg/m . The drag polar C (C ) is shown 

in Fig. 2B and is represented by Eq. (3b), with Kg = 1.235, 

= -2.379, - -5.473. The lift coefficient is subject to 

Ineq. (4), with C = -0.47 and C = -0.21. The maximum lift-to- 

LL LU 

drag ratio modulus is E =0.42 and occurs for C T = -0.47. 

^ max L 

In the numerical experiments, the possibility of a smaller 
lift range is considered via three alternative values for the 
lift coefficient lower bound: 

(i) C LL = -0.47, (ii) C LL = -0.38, (iii) C LL = -0.27. (15) 

For these cases, the values of the maximum lift-to-drag ratio 
modulus are 

(1) E max " °- 42 ' (ii) E max = °' 28 ' (ill) E max " °- 13 (16) 

and occur for C = C . 

L LL 

It must be noted that, for Spacecraft S2, the lift coefficient 

C T = -0.47 corresponds to the angle of attack a = 27 deg; 

C T = -0.38 corresponds to a = 17 deg; C = -0.27 corresponds 
L L 

to a = 10 deg; and C = -0.21 corresponds to a = 7 deg (see Ref. 11). 

L 
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5 . Ideal Optimal Trajectory , Lift Coefficient Unbounded 

An ideal optimal trajectory can be determined analytically 
if one minimizes the total characteristic velocity (10) , subject 
to all of the constraints of Section 2, except the control constraint 
(4) . This can be done via the four-step procedure described below. 

Step 1 . Ignore not only the control constraint (4) , but 
also the differential equations (1) and the static boundary 
conditions (8a) and (9a) . Then, consider the problem of minimizing 
the total characteristic velocity (10) with respect to the terminal 
state components V Q , y Q , , y^ which satisfy the dynamic boundary 
conditions (8b) and (9b) . Clearly, this is a mathematical programming 
problem in which the initial state components V^, y^ and the 
final state components V^, y^ appear separately in both the 
objective function and the constraining relations. Therefore, the 
solution can be obtained by solving two separate mathematical 
programming problems: (i) the minimization of the HEO characteristic 

velocity (10b) with respect to the initial state components V Q , y Q 
which satisfy the constraint (8b); and (ii) the minimization of 
the LEO characteristic velocity (10c) with respect to the final 
state components , y 1 which satisfy the constraint (9b). This 
yields the inertial solutions 


V(> - v »' /[2r 00 /(r 00 + r a n - 

(17a) 

Vi = V./[2r u /(r 11 + r a H, 

(17b) 


and 
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T 0 = 0, 


(18a) 


72 = 0 . 


(18b) 


In turn, because of the inverse transformation relations (7) , the 
inertial solutions (17) -(18) imply the fixed-Earth solutions 


V 0 V * y/ ^ 2r 00 // ^ r 00 + r a^ wr a' 

V ! = V*/[2r 11 /(r 11 + r ft ) ] - u>r a , 


(19a) 

(19b) 


and 


Y 0 = 0, 


72 = 0 . 


(20a) 


(20b) 


Step 2 . Because of the result (20) , it is natural to postulate 
the following altitude/path inclination solutions 


h (t) = h , 

a 


7 ( t) = 0, 


0 < t < 1, 


0 < t < 1, 


(21a) 


(21b) 


which are clearly consistent with the differential equation (la) 
plus the static boundary conditions (8a) and (9a) . 

Step 3. After observing that, because of (21b) , 


7 (t) = 0, 


0 < t < 1, 


( 22 ) 


we see that the differential equation (lc) reduces to 


L = -m(V /r - g ) - m(2wV) , 

cL cl 


( 23 ) 
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which yields the required lift coefficient 

C L = -C 1 [l - (V*/V) 2 ] - C 2 (vyv), (24a) 

where the dimensionless constants , C 2 are given by 


C 1 = 2m/p a Sr a 

11 

CM 

U 

4mw/ p^SV* . 

(24b) 

In particular, the entry and 
coefficient are given by the 

exit values of the required lift 
relations 


C L0 = - C 1 [1 - 

(V./V 0 ) 2 J 

- c 2 <v„/v 0 ). 

(25a) 

C L1 - " C l [1 - 

(V,/ Vl ) 2 ] 

- <V V * / V' 

(25b) 


which must be employed in conjunction with Eqs. (19) . 

Step 4. It remains to verify whether the differential equation 
(lb) can be satisfied. Simple manipulations, omitted for the 
sake of brevity, show that this is precisely the case if the 
velocity-dimensionless time relation V(t) is determined through 
the inverse relation t(V) represented by 


t 


(C 3 A) 



(V*/C D V 2 )dV. 


(26a) 


Here, the dimensional constant (sec) and the flight time x (sec) 
are given by 


C 3 = 2m/p a SV*, 

(26b) 

T = Ca 

r v o 2 

u (v*/c V^dV, 

(26c) 

J J 

V 1 



and the drag coefficient C D is given by Eq. (3b) employed in 
conjunction with Eqs. (24). 
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6 . Properties of the Ideal Optimal Trajectory 

The ideal optimal trajectory determined for lift coefficient 
unbounded is governed by Eqs . (17 ) - ( 26) and satisfies all of the 
constraints of Section 2 except the control inequality constraint 
(4). This trajectory is called grazing trajectory because the 
atmospheric pass is made by flying at constant altitude along the 
edge of the atmosphere. Thus, the complete transfer maneuver 
includes: (i) a descending flight branch connecting an apogee 

located at GEO to a perigee located on the atmospheric edge; 

(ii) a level flight branch flown along the atmospheric edge; and 

(iii) an ascending flight branch connecting a perigee located 
on the atmospheric edge to an apogee located at LEO. 

6.1. Characteristic Velocity . By combining Eqs. (10) and 
(17) -(18), we see that the characteristic velocity of the grazing 
trajectory is given by 

II = AV = AV qo + AV , (27a) 

M 00 = v *^ r a /r 00> - v * (r a /r 00 >/[2r 00 /(r 00 + r a> 1 ' (27b) 

AV i;l = v * / ( r a / r 11 ) - v * (r a /r ll ) ,/[2r ll /(r ll + r a )] * (27c) 

For comparison purposes, consider the Hohmann transfer maneuver, 
which is flown totally in space by navigating in the region bounded 
by r = r^Q and r = r^. For this trajectory, the characteristic 
velocity is given by 

II = AV = AV q0 + AV n , (28a) 


AV 00 V *’ /(r a //r 00 ) V *' / ^ r a / ' r 00 ) ,/ ^ 2r ll //(r 00 + r ll ) 1 ' 


(28b) 
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4 Al “V,/|r a / r n ) + V*/?r a /r u )/l2r 00 /(r 00 + r u >]. (28c) 

Table 1 compares the grazing solution (first line) and the 
Hohmann transfer solution (third line) and shows that the total 
characteristic velocity of the grazing solution is less than 40% 
of that of the Hohmann transfer solution. This illustrates the 
savings in propellant mass which are possible by employing 
aeroassisted orbital transfer techniques instead of Hohmann transfer 
techniques . 

6.2. Peak Altitude Drop . For the grazing trajectory, the 
functional (11) takes the value 

12 = max(h - h) = 0. (29) 

t a 

Therefore, the grazing trajectory simultaneously minimizes the 
peak value of the altitude drop. 

6.3. Peak Dynamic Pressure . For the grazing trajectory, the 

functional (12) takes the value 

13 = max(DP) = (l/2)p V 2 , ( 3 °) 

t cl U * ' 

meaning that the peak dynamic pressure occurs at atmospheric 
entry (t = 0) . Note that the performance index (12) is dominated 
by the density rather than the velocity in the following sense: 
during the atmospheric penetration, the density variations are 
much larger than the velocity variations. Hence, it is felt 
plausible that the grazing trajectory nearly minimizes 
the functional (12) . 

6.4. Peak Heating Rate . For the grazing trajectory, the 
functional (13) takes the value 
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14 = max (HR) = C/ ( P a /p J ( V Q /VJ 3 * 07 , f 31) 

meaning that the peak heating rate occurs at atmospheric entry 
(t = 0) . Once more, because the performance index (13) is 
dominated by the density rather than the velocity, it is felt 
plausible that the grazing traj ectory nearly minimizes the 
functional (13) . 

6.5. Implications . Associated with the grazing trajectory 
(C unbounded) are the following engineering implications for 
a real optimal 'traj ectory (C L bounded). 

(i) The characteristic velocity of the grazing trajectory 
(27) constitutes a lower bound for the characteristic velocity of 
a real optimal trajectory. Therefore, it constitutes a benchmark 
that one should strive to approach with a real optimal trajectory. 

(ii) Because the grazing trajectory simultaneously minimizes 
or nearly minimizes the performance indexes (H) through (14) , 
one surmises that the same general behavior might hold for a 

real optimal trajectory. 

(iii) Along a grazing trajectory, the lift, the centrifugal 
force due to the Earth's curvature, the weight, and the Coriolis 
force due to the Earth's rotation are in static equilibrium. 

(iv) For a grazing trajectory, the velocity decreases 
monotonically between the entry value (19a) and the exit value 
(19b) , and the required lift coefficient increases monotonically 
between the entry value (25a) and the exit value (25b) . For the 
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two spacecraft under consideration, the terminal values of C L 
are given by 


(Si) 

C L0 = - 1921 ‘ 3 ' 

C L1 

( S 2 ) 

C L0 = -750 - 6 ' 

C L1 


-56.6, 


- 22 . 1 , 


(32a) 

(32b) 


with the following meaning: not only is the required lift 
coefficient negative everywhere, but it undershoots by a 
considerable amount the lower bound C . This means that the 

LLi 

grazing trajectory is only an ideal trajectory, hence not flyable. 

(v) To sum up, while the grazing trajectory was computed neglecting 

the two-sided inequality constraint (4) , the nature of the grazing 

solution is such that the upper bound inequality C < C is 

satisfied everywhere, while the lower bound inequality C > C 

Li L.L 

is violated everywhere. Even though the grazing trajectory is 
not flyable, this suggests the idea that, for a real optimal 
trajectory, the control distribution might be of the form C L = C LL . 

This idea is explored in the following section. 
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7 . Real Optimal Tra j ectory , Lif t Coefficient Bounded 

A real optimal trajectory can be determined numerically by 
accounting for all of the constraints of Section 2, including 
the control constraint (4) . Two alternative procedures are now 
presented . 

7.1. Solution via TPBVP . The grazing trajectory of Sections 5-6 
satisfies all of the constraints of Section 2, except the two-sided 
inequality constraint (4) . Specifically, the lower bound inequality 
C T > C TT is violated everywhere. This suggests the idea that the 
control distribution of the real optimal trajectory might be of 
the form C L = . 

If the hypothesis C T = C TT is made, the resulting trajectory 

Li XjLi 

is described by the differential system 


h = t [V siny] , 

(33a) 

• 

V = x[-D(h,V,C )/m - gsiny] , 

(33b) 

y = x [ L(h,V,C )/mV + (V/r - g/V)cosy + 2u] , 
LL 

(33c) 

T = 0, 

(33d) 


which must be solved in conjunction with the initial conditions 


h 0 = V 


(34a) 


<2r 00 - 2r 00 r a )V * 2 - r 00 (V 0 2 + 2wr a V 0 COSY 0 + “M’ 


+ r a 2 (V 0 coSY 0 + ojr a ) 2 = 0, 


(34b) 


and the final conditions 


\ 



21 


AAR-249 


h, = h , (35a) 

L cl 

(2r ll ‘ 2r U r a )V . 2 - r ll (V l + 2t0r a V l cos A + “M’ 

2 2 

+ r (V.cosy. + tor ) =0. (35b) 

a 1 '1 a 

The following comments are pertinent: 

(i) the system (33) is a formal modification of the system 
(1) after accounting for the hypothesis C T = C and the fact 
that, since x is a parameter, it can be represented via the 
differential equation x = 0; this is done only for explanation purposes 

(ii) the boundary conditions (34b) and (35b) are a formal 
modification of the boundary conditions (8b) and (9b) after accounting 
for the transformation relations (5)— (6) ; 

(iii) the problem represented by Eqs. (33) -(35) is a two- 
point boundary -value problem (TPBVP) in which the unknowns are 

the functions h(t), V(t), y (t) , x (t) ; since there are four differential 
equations and four boundary conditions, one surmises that a 
solution might exist; however, because the TPBVP (33) -(35) is 
nonlinear, the existence of a solution must be confirmed by numerical 
tests . 

7.2. Solution via SGRA . An alternative to solving the previous 
TPBVP is to employ the sequential gradient-restoration algorithm 
(SGRA, Refs. 16-18) in order to minimize the functional (10) with 
respect to the state variables h(t), V(t), y (t) , the control 
variable C (t) , and the parameter x which satisfy the differential 

Li 

system (1)— (3) , the control constraint (4), and the boundary 
conditions (8) -(9) rewritten in the form (34) -(35) after replacement 
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of inertial quantities with Earth -fixed quantities. Note that 
SGRA is an iterative algorithm which produces a sequence of 
feasible solutions, each characterized by a lower value of the 
functional (10) . While SGRA is available in both primal form 
and dual form (Refs. 17-18), the primal form is more suitable 
to the solution of hypervelocity flight problems. Hence, it is 
employed here. 

7.3. Numerical Results . In computing numerical solutions 
via SGRA, the experimental data of Section 4 were employed. Hence, 
a GEO-to-LEO optimal aeroassisted orbital transfer was determined 
for both Spacecraft SI (see Refs. 13-14) and Spacecraft S2 (see 
Ref. 11) . For each spacecraft, three different values were considered 
for the lift coefficient lower bound, specifically, the values 
(14) for Spacecraft SI and the values (15) for Spacecraft S2. 

The numerical results are shown in Tables 1-2 and Figs. 3-4. 

Table 1 compares the characteristic velocity of the real 
optimal trajectory (second line) with that of the grazing trajectory 
(first line) and the Hohmann transfer trajectory (third line) . 
Clearly, the characteristic velocity of the real optimal trajectory 
is within 1% of that of the grazing trajectory and is considerably 
below (by almost 60%) that of the Hohmann transfer trajectory. 

Table 2 presents summary results for the real optimal trajectory 
of Spacecraft SI and Spacecraft S2. In Table 2, the following 
quantities are shown: the entry values of the state variables; the 
exit values of the state variables; the minimum altitude, the 
peak dynamic pressure, the peak heating rate, and the flight time; 


\ 
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the components of the characteristic velocity and the total 
characteristic velocity. As C T increases (namely, as the lift 

J-jXj 

coefficient range C - C decreases), the real optimal trajectory 
is characterized by a steeper entry, lower minimum altitude, higher 
peak dynamic pressure, and higher peak heating rate. However, 
the total characteristic velocity and its components are almost 
independent of C , hence almost independent of the lift 
coefficient range. 

Figures 3-4 refer to the real optimal trajectories of Spacecraft 

SI and Spacecraft S2. The following quantities are shown as 

functions of the dimensional time T: the altitude h; the velocity 

V; the path inclination y; the lift coefficient C ; and the density 

p. The results confirm the idea that, for the real optimal trajectory, 

the lift coefficient is constant and equal to the lower bound. 

As C T _ increases (namely, as the lift coefficient range C T r7 - C 

LL LU LL 

decreases), the real optimal trajectory changes; in particular, the 
minimum altitude changes. However, while the changes in minimum 
altitude are relatively small, they translate into larger density 
changes, and hence larger changes in max(DP) and max(HR). 


1 
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8 • Properties of the Real Optimal Trajectory 

The real optimal trajectory determined for lift coefficient 
bounded satisfies all of the constraints of Section 2, including 
the control inequality constraint (4). This trajectory is called 
nearly-grazing trajectory, because the atmospheric pass is 
characterized by atmospheric penetration with the smallest possible 
entry angle, followed by flight at the lift coefficient lower 
bound. Thus, the complete transfer maneuver includes: (i) a 

descending flight branch connecting an apogee located at GEO to 
the atmospheric entry point; (ii) an atmospheric flight branch 
including a shallow descent followed by a shallower ascent; and 
(iii) an ascending flight branch connecting the atmospheric exit 
point to an apogee located at LEO. 

8.1. Characteristic Velocity . If we combine Eqs. (8) -(10) 
and eliminate the terminal velocities, we see that the characteristic 
velocity components of the nearly-grazing trajectory can be expressed 
in terms of the terminal path inclinations as follows: 

^00 = v * / ' r a /r 00> - v * (r a /r OO )/T2r OO (r OO * r a ,cos2 V 

<r oo - <36a) 

= V,^r a /r u ) - V.tr a /r 11 )/t2r 11 (r 11 - r^cos 2 ^/ 

(r^ - r^cos^y, ) ] . (36b) 

11 ci 1. 

By comparison, we observe that, for the grazing trajectory 



\ 
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AV 


OOG 


AV 


11G 


V*/Tr a /r 

V*^r a /r 


00 ' 


lr 


V * (r a /r 00 )/f2r 00 /(r 00 + r a n ' 

(37a) 

V * (r a /r ll )/t2r ll /(r ll + r a )1 ' 

(37b) 


Let the right-hand sides of (36) be expanded in a Taylor series/ 
and let the following approximations be employed: 


Y 0 = Y 0 , (38a) 

h = Y 1 . (38b) 


Then, the following relations can be established between (36) and 
(37): 


4V 00 

av oog + A 0 Y 0 ' 

(39a) 

fiVu 

A ^11G + A 1 Y 1 ' 

(39b) 


with 

A 0 = <V,//2)[r a /(r 00 - r a )][r 00 /(r 00 + r a )] 3/2 , 

\ - (V,/A2)[r a /(r 11 - r a IMr u /(r u + r a )] 3/2 . 

For the GEO-to-LEO transfer of Section 4, the constants A^ , A^ 
take the values 


A q = 0.8138 km/sec, (41a) 

A 1 = 62.03 km/sec. (41b) 

Equations (39) show that the characteristic velocity components of 
a nearly-grazing trajectory can be contained by making the entry 
and exit path inclinations as small as possible. For tangential entry 
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and exit (y Q = y^ = 0) , the characteristic velocity components 
of the nearly-grazing trajectory reduce to those of the grazing 
trajectory. For nontangential entry and exit (y^ ? 0 , y^ / 0), 
the characteristic velocity components of the nearly-grazing 
trajectory are always larger than those of the grazing trajectory. 
However, the relative differences are small, of order less 
than 1% for the examples of Tables 1-2. 

3.2. Terminal Velocities . Equations (8b) and (9b) can be 
solved to obtain the terminal velocities of the nearly-grazing 
trajectory in terms of the terminal path inclinations as follows: 

= v . /i2r 00 (r 00 ■ r a ,/(r 00 • r a 2 cos 2 V 1 ' (42a> 

= V vt /(2r 11 (r 11 - r a )/(i'i 1 ~ r a 2 cos 2 y 1 )]. (42b) 

By comparison, we observe that, for the grazing trajectory 

6 0 = = o), 


0G V *’ / ~ [2r 00 /(r 00 

+ r )] , 

d 

(43a) 

1G = v V 2r n/( r u 

+ r ) ] . 

d 

(43b) 


Let the right-hand sides of (42) be expanded in a Taylor series, 
and let the approximations (38) be employed. Then, the following 
relations can be established between (42) and (43) ; 


V 0 V 0G VO ' 


V 1 = V 1G - Vl ' 


(44a) 


\ 


(44b) 
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which approximately imply that 
V 0 = V 0G - Vo"' 


V 1 = V 1G - Vl ' 


(45a) 

(45b) 


with 


2 // 2 2 
'a. r 00 - r a 


)]/[r 00 / (r 00 

+ r )] , 

d 

(46a) 

) ] (r ll 

+ r ) ] . 

d 

(46b) 


For the GEO-to-LEO transfer of Section 4, the constants B^, 
take the values: 


0.1254 

km/sec, 

(47a) 

60.09 

km./ sec . 

(47b) 


Equations (44) -(45) show that the terminal velocities of a nearly- 
grazing trajectory can be made close to the terminal velocities 
of a grazing trajectory by making the entry and exit path 
inclinations as small as possible. For tangential entry and exit 
( Yq = y = 0) , the terminal velocities of the nearly-grazing 
trajectory reduce to those of a grazing trajectory. For nontangential 
entry and exit (y q ? 0, f 0) , the terminal velocities of the 
nearly-grazing trajectory are always smaller than those of the 
grazing trajectory. However, the relative differences are small, 
of order 10 -4 (entry velocities) or order 10 3 (exit velocities) 
for the examples of Tables 1-2. 


\ 
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8.3. Perigee Altitude . Let us consider the trajectory 
described by the spacecraft in the absence of aerodynamic forces 
and evaluate the energy constant twice, once in terms of the entry 
conditions and once in terms of the exit conditions. For the 
nearly-grazing trajectory, this yields the relations 


E 0 = V 0 2/2 ‘ p/r a = ■‘ ,/(r 00 + r P0>' (48a) 

E x = V 1 2 /2 - y/r a = -ii/(r i;L + r pl ) , (48b) 


in which r . r , are the perigee radii associated with the 
PO Pi 

entry and exit conditions, respectively. The corresponding relations 
for the grazing trajectory are 


E 0G ^G^ 2 y//r a -y// ^ r 00 + r a^ ' 


E 1G = V 1G /2 - v/r a “ _li/(r n + r a>’ 


(49a) 

(49b) 


Upon taking differences, we see that 


<*0 

ii 

CM 

\ 

*"b 

CM O 
?> 

1 

-u (r 

a 

r p0^/( r 00 

+ r P0 )(r 00 + r a > ' 

(50a) 

<°1 

- 0 1G> /2 * 

-!i< r a 

- r P1 >/<r u 

+ r Pl> (r ll + r a’ ' 

(50b) 

we 

observe that 




r 

a 

■ r P0 = h a 

h P0' 



(51a) 

r 

a 

- r Pl " h a 

- h pl . 



(51b) 


and employ the approximations 


<r 00 + r P0> (r 00 + r a> £ (r 00 + r a> ' 


(52a) 


\ 


(r ll + r Pl J (r ll + r a ) = (r xl + r a ) 2 . 


(52b) 
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With this understanding, Eqs . (50) can be rewritten as 


< V 0 V 0G^ 2 y ^ h a h P0^ r 00 + r a^ ' 

(V 2 - V 2 g )/2 = -U(h a - h pl )/(r u + r a ) 2 . 


(53a) 

(53b) 


Because y = r V , Eqs. (53) imply that 

3 . * 


h a h P0 " (V 0 (r 00 + r a } //2r a V * ' 

h a " h Pl = " ( ^1 - ^1G> (r ll + r a )2/2r a V * 2 ' 


(54a) 

(54b) 


Note that, in the light of (44), 


^0 " V 0G 2B 0^0G Y 0 ' 


*1 - *Ig =" 2b iVi 2 ' 


(55a) 

(55b) 


Hence, the altitude drops of the ideal perigees can be written as 


h a h P0 C 0 Y 0 ' 


h a h Pl C 1 Y 1 ' 


(56a) 

(56b) 


where 


C 0 = r a r 00 /<r 00 - r a>' 


C 1 = r a r ll / (r ll ’ r a> 


(57a) 

(57b) 


For the GEO-to-LEO transfer of Section 4, the constants C Q , C ^ 
take the values 


Cq = 7682 km, 


(58a) 


4 



30 


AAR-249 


C ± = 207565 km. (58b) 

Equations (56) show that the altitude drop of the ideal perigee 
of a nearly-grazing trajectory can be made small by making the 
entry and exit path inclinations as small as possible. For 
tangential entry and exit (y^ = y^ = 0) , the altitude drop of the 
ideal perigee of a nearly-grazing trajectory reduces to that of a 
grazing trajectory, which is zero. For nontangential entry and 
exit (y Q ^ 0, y^/0), the altitude drop of the ideal perigee 
is nonzero. In particular, Eq. (56a) yields a good approximation 
to the minimum altitude of a nearly-grazing trajectory (within 1 km) . 

8.4. Minimum Altitude . The perigee altitude determined with 
either of Eqs. (56) constitutes an approximation to h ^ n , the 
minimum altitude of the nearly-grazing trajectory. The approximation 
is more accurate if Eq. (56a) is used instead of Eq. (56b) , due 
to the fact that the descending branch of the nearly-grazing 
trajectory is shorter in time than the ascending branch. Hence, 
the effects due to the aerodynamic forces are smaller in the 
descending branch than in the ascending branch. 

Alternatively, one might regard the left-hand sides of Eqs. (56) 

as approximations to max (h - h) , the maximum altitude drop of 

t 3 

the nearly-grazing trajectory. To contain maxdi^-h), one must make 
the right-hand sides of Eqs. (56) as small as possible. Hence, 
the entry and exit path inclinations should be as small as 
possible. One surmises that the solution of Problem (P2) , minimum 
peak altitude drop, is nearly the same as the solution of 
Problem (Pi) , minimum energy. 
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8.5. Dynamic Pressure . During the atmospheric pass of an 
AOT trajectory, both the density and the velocity change, and 
hence the dynamic pressure (12b) changes. The ratio of the 
instantaneous dynamic pressure to the entry dynamic pressure 
is given by 

DP/(DP) q = (p/p a ) (V/V Q ) 2 . 

While the velocity ratio V/V^ is of 0(1) , the density ratio p/p 
3 4 

is of 0(10 ) to 0(10 ), depending on the particular trajectory. 
Therefore, to contain the dynamic pressure, one must act primar 
on the density, hence the altitude, hence the drop of altitude. 
One surmises that the solution of Problem (P3) , minimum peak 
dynamic pressure, is nearly the same as the solution of 
Problem (P2) , minimum peak altitude drop; hence, it is nearly 
the same as the solution of Problem (Pi) , minimum energy. 

8.6. Heating Rate . During the atmospheric pass of an AOT 
trajectory, both the density and the velocity change, and hence 
the heating rate (13b) changes. The ratio of the instantaneous 
heating rate to the entry heating rate is given by 

8 07 

HR/ ( HR) 0 = /(p/p a )(V/V Q ) ' . 

While the velocity ratio V/Vq is of 0(1), the density ratio p/p 
3 4 

is of 0(10 ) to 0(10 ), depending on the particular trajectory. 
Therefore, to contain the heating rate, one must act primarily 
on the density, hence the altitude, hence the drop of altitude. 
One surmises that the solution of Problem (P4) , minimum peak 
heating rate, is nearly the same as the solution of , 


(59) 


iy 


(60) 
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Problem (P2) , minimum peak altitude drop; hence, it is 

nearly the same as the solution of Problem (Pi) , minimum energy. 

8.7. Lift Coefficient . For a nearly-grazing trajectory, 
both y and y are small along a major portion of the trajectory; 
in particular, they are small in the region past the minimum 
altitude point. For y and y small, Eq. (1c) can be rewritten as 

L = (l/2)C T pSV^ = -m(V^/r - g) - m(2wV). (61) 

Li 

Because the atmospheric pass is executed at supercircular velocities, 

2 

the centrifugal force due to the Earth's curvature m(V /r) always 

exceeds the weight mg; also, the Coriolis force m(2uV) is directed 

upward. Because both terms on the right-hand side of (61) are 

negative, the required lift is negative; hence, the required lift 

coefficient is negative. To maximize the instantaneous flight 

altitude, the density p should be as small as possible; hence, 

the lift coefficient C T should be as negative as possible. To 

sum up, the control C T = C TT allows the flight altitude to be 

L LL 

as high as possible, and this is why the solutions of Problems (PI) 
through (P4) are nearly the same. 

8.8. Summary . The real optimal trajectory determined in the 

presence of the control constraint (4) is a nearly-grazing 

trajectory and is flown with the lift coefficient at the lower 

bound, C T = C T T . The control distribution C T = C T T not only 
L LL -L- 1 JjJ-j 

minimizes the characteristic velocity of Problem (Pi) , but 
simultaneously nearly minimizes the peak altitude drop of Problem 
(P2) , the peak dynamic pressure of Problem (P3) , and the peak 
heating rate of Problem (P4) . 
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9 . Conclusions 

This paper is concerned with the optimization of 
trajectories for coplanar, aeroassisted orbital transfer (AOT) 
from a geosynchronous Earth orbit (GEO) to a low Earth orbit 
(LEO) . It is assumed that the initial and final orbits are 
circular, that the gravitational field is central and is governed 
by the inverse square law, and that two impulses are employed, 
one at GEO exit and one at LEO entry. During the atmospheric 
pass, the trajectory is controlled via the lift coefficient in 
such a way that the total characteristic velocity is minimized. 

First, an ideal optimal trajectory is determined 
analytically for lift coefficient unbounded. This trajectory is 
called grazing trajectory, because the atmospheric pass is made 
by flying at constant altitude along the edge of the atmosphere 
until the excess velocity is depleted. For the grazing 
trajectory, the lift coefficient varies in such a way that the 
lift, the centrifugal force due to the Earth’s curvature, the 
weight, and the Coriolis force due to the Earth's rotation are in 
static balance. Also, the grazing trajectory minimizes the total 
characteristic velocity and simultaneously nearly minimizes the 
peak values of the altitude drop, the dynamic pressure, and the 
heating rate. 

Next, starting from the grazing trajectory results, a real 
optimal trajectory is determined numerically for lift coefficient 
bounded from both below and above. This trajectory is 
characterized by atmospheric penetration with the smallest 
possible entry angle, followed by flight at the lift coefficient 
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lower bound. It is called nearly-grazing trajectory, in that its 
control distribution is the closest feasible control to that of 
the grazing trajectory. Consistently with the grazing trajectory 
behavior, the nearly-grazing trajectory minimizes the total 
characteristic velocity and simultaneously nearly minimizes the 
peak values of the altitude drop, the dynamic pressure, and the 
heating rate. 



35 


AAR-249 


References 

1. BOSTON, P. J. , Editor, The Case for Mars , Univelt, San Diego, 
California, 1984. 

2. REIBER, D. B., Editor, The NASA Mars Conference , Univelt, 

San Diego, California, 1988. 

3. REA, D. G., CRAIG, M. K. , CUNNINGHAM, G. E., and CONWAY, H. L. , 
An International Mars Exploration Program , Acta Astronautica , 

Vol . 22, Nos. 1-12, pp. 255-267, 1990. 

4. MEASE, K. D., and VINH, N. X., Minimum-Fuel Aeroassisted 
Coplanar Orbit Transfer Using Lift Modulation , Journal of 
Guidance, Control, and Dynamics, Vol. 8, No. 1, pp. 134-141, 

1985. 

5. MIELE, A., and VENKATARAMAN , P., Optimal Trajectories for 
Aeroassisted Orbital Transfer , Acta Astronautica, Vol. 11, 

Nos. 7-8, pp. 423-433, 1984. 

6. MIELE, A., and BASAPUR, V. K. , Approximate Solutions to Minimax 
Optimal Control Problems for Aeroassisted Orbital Transfer , 

Acta Astronautica, Vol. 12, No. 10, pp. 809-818, 1985. 

7. MIELE, A., BASAPUR, V. K., and MEASE, K. D. , Nearly-Grazing 
Optimal Trajectories for Aeroassisted Orbital Transfer , 

Journal of the Astronautical Sciences, Vol. 34, No. 1, pp. 3-18, 

1986. 

8. MIELE, A., BASAPUR, V. K., and LEE, W. Y. , Optimal Trajectories 
for Aeroassisted, Coplanar Orbital Transfer , Journal of 
Optimization Theory and Applications, Vol. 52, No. 1, 

pp. 1-24, 1987. 



36 


AAR-249 


9. VINH, N. X., JOHANNESEN, J. R. , MEASE, K. D. , and HANSON, J. M. , 
Explicit Guidance of Drag-Modulated Aeroassisted Transfer 
between Elliptical Orbits , Journal of Guidance, Control, and 
Dynamics, Vol. 9, No. 2, pp. 274-280, 1986. 

10. LEE, B. S., and GRANTHAM, W. J., Aeroassisted Orbital 
Maneuvering Using Lyapunov Optimal Feedback Control , Journal 
of Guidance, Control, and Dynamics, Vol. 12, No. 2, pp. 237- 
242, 1989. 

11. ANONYMOUS, N. N. , Aeroassisted Flight Experiment : 

Preliminary Design Document , NASA Marshall Space Flight 
Center, 1986. 

12. GAMBLE, J. D., CERIMELE , C. J. , MOORE, T. E. , and HIGGINS, J. , 
Atmospheric Guidance Concepts for an Aeroassist Flight 
Experiment , Journal of the Astronautical Sciences, Vol. 36, 

Nos. 1-2, pp. 45-71, 1988. 

13. MIELE, A., WANG, T. , and LEE, W. Y. , Optimization and Guidance 
of Trajectories for Coplanar, Aeroassisted Orbital Transfer , 

Rice University, Aero-Astronautics Report No. 241, 1989. 

14. MIELE, A., and WANG, T. , Gamma Guidance of Trajectories for 
Coplanar, Aeroassisted Orbital Transfer , Rice University, 
Aero-Astronautics Report No. 246, 1990. 

15. NOAA, NASA, and USAF, US Standard Atmosphere, 1976, US 
Government Printing Office, Washington, DC, 1976. 

16. MIELE, A., WANG, T., and BASAPUR, V. K. , Primal and Dual 
Formulations of Sequential Gradient-Restoration Algorithms 
for Trajectory Optimization Problems , Acta Astronautica , 

Vol. 13, No. 8, pp. 491-505, 1986. 


\ 



37 


AAR-249 


17. MIELE, A., and WANG, T. , Primal-Dual Properties of Sequential 
Gradient-Restoration Algorithms for Optimal Control Problems , 
Part 1, Basic Problem , Integral Methods in Science and 
Engineering, Edited by F. R. Payne et al. Hemisphere 
Publishing Corporation, Washington, DC, pp. 577-607, 1986. 

18. MIELE, A., and WANG, T., Primal-Dual Properties of Sequential 
Gradient-Restoration Algorithms for Optimal Control Problems , 
Part 2, General Problem , Journal of Mathematical Analysis 
and Applications, Vol . 119, Nos. 1-2, pp. 21-54, 1986. 



38 


AAR-249 


List of Tables 

Table 1A. Comparison of characteristic velocities (km/sec) 
Spacecraft SI. 

Table IB. Comparison of characteristic velocities (km/sec) 
Spacecraft S2. 

Table 2A. Optimal trajectory results. Spacecraft SI. 

Optimal trajectory results. Spacecraft S2. 


Table 2B. 



39 


AAR-249 


Table 1A. Comparison of characteristic velocities (km/sec) , 
Spacecraft SI. 


— 

Trajectory 

Characteristic 

velocity 

C=-0.90 

LiJj 

C =-0 . 70 

i-ilJ 

C _ =-0 . 50 

J— i JLi 


Grazing 

A ^oo 

1.4857 

1.4857 

1.4857 

— 

Optimal 

^00 

1.4899 

1.4901 

1.4903 


Hohmann 

o 

o 

?> 

o 

1.4638 

1.4638 

1.4638 


Grazing 

A Vn 

0.0615 

0.0615 

0.0615 


Optimal 

AV U 

0.0673 

0.0675 

0.0678 

— 

Hohmann 

AV 11 

2.4173 

2.4173 

2.4173 


Grazing 

AV 

1.5472 

1.5472 

1.5472 


Optimal 

AV 

1.5572 

1.5576 

1.5581 


Hohmann 

AV 

3.8810 

3.8810 

3.8810 
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Table IB. Comparison of characteristic velocities (km/sec), 


Spacecraft S2. 

Trajectory 

Characteristic 

r- 

• 

o 

i 

II 

A 

h-q 

U 

C =-0.38 

Lb 

c liT-°- 27 


velocity 




Grazing 

A \o 

1.4857 

1.4857 

1.4857 

Optimal 

^00 

1.4896 

1.4897 

1.4899 

Hohmann 

^00 

1.4638 

1.4638 

1.4638 

Grazing 

A ^ll 

0.0615 

0.0615 

0.0615 

Optimal 

AVn 

0.0684 

0.0702 

0.0733 

Hohmann 

AVu 

2.4173 

2.4173 

2.4173 

Grazing 

AV 

1.5472 

1.5472 

1.5472 

Optimal 

AV 

1.5580 

1.5599 

1.5632 

Hohmann 

AV 

3.8810 

3.8810 

3.8810 
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Table 

2A. Optimal 

trajectory 

results. Spacecraft SI. 

Quantity 

C =-0. 90 

LiLj 

C_ =-0.70 
LiXj 

C =-0 . 50 
LLi 

Units 

h 0 

120.00 

120.00 

120.00 

km 

V 0 

9.8371 

9.8371 

9.8372 

km/ sec 

v 0 

10.3097 

10.3096 

10.3096 

km/sec 

Y o 

-4.297 

-4.388 

-4.508 

deg 

^0 

-4.100 

-4.187 

-4.301 

deg 

h l 

120.00 

120.00 

120.00 

km 

V 1 

7.4146 

7.4145 

7.4141 

km/ sec 

*1 

7.8884 

7.8883 

7.8880 

km/sec 

Y 1 

0.592 

0.601 

0.616 

deg 

h 

0.557 

,0.565 

0.579 

deg 

min (h) 

75.35 

73.46 

70.96 

km 

max (DP) 

1675 

2245 

3264 

N/m^ 

max (HR) 

177.1 

206.3 

251.2 

W/cm^ 

T 

1872.2 

2172.1 

2547.6 

sec 

o 
o 
? > 
<1 

1.4899 

1.4901 

1.4903 

km/ sec 

iVn 

0.0673 

0.0675 

0.0678 

km/sec 

AV 

1.5572 

1.5576 

1.5581 

km/sec 
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Table 

2B. Optimal 

trajectory 

results. Spacecraft 

. S2. 

Quantity 

C LL =-°- 47 

C =-0.38 

Li Li 

C LL =_0 ’ 27 

Units 

h 0 

120.00 

120.00 

120.00 

km 

V 0 

9.8371 

9.8371 

9.8371 

km/ sec 

v 0 

10.3097 

10.3097 

10.3097 

km/ sec 

Y o 

-4.167 

-4.215 

-4.293 

deg 


-3.976 

-4.022 

-4.096 

deg 

h l 

120.00 

120.00 

120.00 

km 

V 1 

7.4136 

7.4119 

7.4089 

km/ sec 

*1 

7.8874 

7.8857 

7.8827 

km/ sec 


0.644 

0.721 

0.841 

deg 

*1 

0.605 

0.678 

0.790 

deg 

min (h) 

78.52 

78.09 

77.25 

km 

max (DP) 

986 

1034 

1153 

N/m 2 

max (HR) 

131.6 

132.7 

138.2 

W/cm 2 

T 

1069.1 

856.4 

691.2 

sec 

D> 

O 

O 

1.4896 

1.4897 

1.4899 

km/ sec 

AVn 

0.0684 

0.0702 

0.0733 

km/sec 

AV 

1.5580 

1.5599 

1.5632 

km/ sec 
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